Atomic structure and electronic properties of the cleavage InAs(llO) surface 
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The electronic properties of the cleavage InAs(llO) surface are studied using a semi-empirical 
tight-binding method which employs an extended atomic- like basis set. We calculate the surface 
electronic band structure which is studied as a function of the structural parameters of the surface. 
We describe and discuss the electronic character of the surface electronic states and we compare 
with other theoretical approaches, and with experimental observations. Finally, we discuss which 
atomic model better resembles the available experimental data. 



I. INTRODUCTION 

Low-dimensional quantum structures have shown to 
have unique optical and electronic properties, being im- 
portant in the fabrication of new opto-electronic devicesi. 
III-V semiconductor compounds have demonstrated to 
be important in the development of new light emitters 
and detectors. In particular, Indium compounds, like 
InAs and InSb, are good infrared detectors due to their 
small bulk-band gap. Most of modern opto-electronic 
devices are fabricated growing monoatomic layers on a 
semiconductor substrate using epitaxial techniques. The 
epitaxial growth depends on the ambient conditions as 
well as in the atomic and electronic properties of the sur- 
face of the substrate. Therefore, the continuous develop- 
ment of semiconductor heterostructures and other struc- 
tures at the nanometer scale, makes evident the need for 
a detailed theoretical and experimental understanding of 
semiconductor surfaces. 

A lot of theoretical work has been done to determine 
the electronic properties of Gallium compounds, and 
their surfaces^. On the other hand, only a few attempts 
have been made to characterize the InAs surfaces. The 
(110) surface is the natural cleavage of zincblende crys- 
tals, and it is a non-polar surface which contains equal 
number of cations and anions in its unit cell, showing a 
partly ionic bonding. In this work, we are interested in 
the InAs(llO) surface. 

Most of the theoretica&i^iSii studies about InAs(llO) 
do not provide a way for directly compare with exper- 
imenta^iQi^°i"i^^i^3i^'^i^-'^i^^i^^i^^ Furthermore, the avail- 
able experimental measurements are not enough to com- 
pletely elucidate the atomic structure and electronic 
properties of InAs(llO). For example, Andersson and 
collaborators^ found the energies of occupied-surface 
states at high-symmetry points using photoemission tech- 
niques. Six years later, the same group measured 
again the occupied-surface states at the same high- 
symmetry pointsS, founding a systematic shift of energy 
of —0.15 eV with their previous measurcmcntE*. Inde- 
pendently, Swantson et al}^ also measured the energies 
of occupied-surface states at high-symmetry points and 
they found differences up to 0.5 eV with those reported 



by Andersson and collaborators*'^. In general, the inter- 
pretation of photoemission spectra has been difficult to 
do because of the small bulk-band gap of InAs. 

Theoretically, the electronic structure and atomic po- 
sitions of InAs(llO) were calculated using a quantum- 
molecular dynamics^ based in a semi-empirical tight- 
binding (TB) approach. Almost a decade latter, an 
ah initio quantum-molecular dynamics^- was performed. 
The reported atomic structure and electronic surface 
states differ between TB and ah initio calculations, and 
also differ with the available experimental measurements. 
This fact is because the semi-empirical calculations were 
performed using an atomic reconstruction that was not 
fully relaxed, while ah initio calculations were performed 
using a DFT-LDA with a plane-wave basis set that de- 
pend on the cut-off energy to get accurate surface states, 
and present systematic errors to determine the precise 
electronic energy levels. In summary, there is no a con- 
sensus on the physical properties of the InAs(llO) sur- 
face. 

In this work, we present the electronic structure of 
InAs(llO) employing a semi-empirical TB formalismii2iS£ 
and using the atomic coordinates obtained from ah ini- 
tio quantum-molecular dynamics^. The use of the fully 
relaxed slab coordinates guarantees that the calculated 
electronic properties include all the subtle effects of 
surface-induced strain and appropriate geometry. The 
TB approach allows us to analyze in detail the electronic 
structure and compare our calculations with available ex- 
perimental data. Our calculated surface electronic band 
structure is very similar to that calculated by Shkrebtii 
and collaborators^, however, in that work they focused to 
study the electronic properties of the adsorption of anti- 
mony on III-V(llO) surfaces, and only a rough analysis of 
the clean InAs(llO) surface was done, and a comparison 
with experiments and theoretical results was not done. 



II. THEORETICAL METHOD 

The III-V(llO) semiconductor surfaces relax in such 
a way that the surface cation atom moves inwards the 
surface into an approximately planar configuration, with 
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a threefold coordination with its first-neighbors anion 
atoms. The topmost anion atom moves outward to the 
surface, showing a pyramidal configuration with its three 
first-neighbors cation atoms2*2i. The geometric parame- 
ters that describe the relaxation of the surface atoms of 
III-V(llO) semiconductor surfaces, scale linearly with the 
bulk lattice constant^^. In particular, Alves et al. found^ 
that for the InAs(llO) surface the pyramidal angle at the 
anion, labeled by a, is ~ 90°, the in-plane angle [3 has 
values close to the tetrahedral bond angle ~ 109.47°, and 
the planar angle at the cation, labeled by 7, is ~ 120°. 
For the ideal surface the values for a, (3 and 7 correspond 
to those angles of tetrahedral bonds, 109.47°. 

The relaxed InAs(llO) surface is shown in Fig. ^ In 
Fig. n (a) we show the top view of a surface unit cell 
that contains one In atom (cation), and one As atom 
(anion) per atomic layer. The open circles correspond 
to As atoms while black circles show In atoms. The pa- 
rameter Co is the bulk lattice constant and do = ao/2\/2. 
In Fig. n (b) we show a side view with only the three 
outermost atomic layers of the slab. Here, we define the 
structural parameters associated to the surface relaxation 
whose values are given in Table |I| In Fig. ^ (c) we show 
the corresponding Two-Dimensional Irreducible Brillouin 
Zone (2DIBZ). 

In our calculations, the non-polar InAs (110) surface 
was modeled using a slab of 50 atoms, yielding a free re- 
constructed surface on each face of the slab. The thick- 
ness of the slab is large enough to decouple the surface 
states at the top and bottom surfaces of the slab. Peri- 
odic boundary conditions were employed parallel to the 
surface of the slab to effectively model an infinite two- 
dimensional crystal system. The atomic coordinates were 
taken from Ref. 3, and are given in Table ID We have 
performed calculations with all the structural parame- 
ters in TableQl however, those corresponding to the Den- 
sity Functional Theorji^ (DFT) with an energy cutoff of 
18 Ry are the ones that best resemble the available exper- 
imental datagi'^ii°i"ii^ii3ii4,i5,i6,i7,i8^ YVe calculate the 

electronic level structure of the slab using a well known 
parameterized TB approach with a sp'^s* orbital-like ba- 
sis, within a first-neighbor interaction approaches. This 
wave function basis provides a good description of the 
valence and conduction bands of cubic semiconductors. 
This TB approximation has been applied to calculate 
the electronic and optical properties of a variety of semi- 
conductor surfaces, including other III-V compounds?^. 
The TB parameters are taken to be the same as those 
of VoglSfi for the bulk but they are scaled by a factor 
of (D/d)^, where d is the bond length of any two first- 
neighbor atoms, and D = -y3ao/4i^ These changes to 
the original bulk parameters provide an excellent descrip- 
tion of the electronic structure, as compared to experi- 
mental measurements. 

The atomic structure of the surface region is intimately 
related to its electronic structure. Experimentally, the 
electronic structure can be determined by means of elec- 
tron spectroscopies like photoemission (PE), inverse pho- 



toemission (IPE), and Scanning Tunneling Spectroscopy 
(STS). These techniques are sensitive to the surface's fea- 
tures and electronic properties due to reconstructions or 
adsorption events. We present in Section III the surface 
electronic band structure and the local density of elec- 
tronic states of the reconstructed or relaxed InAs(llO) 
surface. 



III. RESULTS 

We show the surface electronic band structure along 
high-symmetry points of the 2DIBZ of the reconstructed 
InAs(llO) surface in Fig. |2 The projected bulk elec- 
tronic states are shown in tiny black dots, while the sur- 
face electronic states are shown in large black dots. The 
calculated Fermi energy level, Ep, is at 1.1 eV above the 
Valence-Band-Maximum (VBM), while the Minimum of 
the Conduction Band (MBC) is E = 0.6 eV from 
the VBM at the high-symmetry point P. We denote the 
surface electronic states using the labels Ai and Ci asso- 
ciated to the surface anions and cations, respectively, as 
introduced by Chelikowsky and Cohen24. 

Below the VBM we found four well-defined occupied 
surface electronic states denoted by A5, ^3, A2 and C2. 
The Ac, surface states correspond to the dangling bonds 
of the As atoms located at the first atomic layer. The 
states form a band from the high-symmetry point X to 
the point X', going through the high-symmetry point M 
in the 2DIBZ. This band has a minimum at X with an 
energy of -1.20 eV and disperse upwards towards the F 
point. From X, the band also disperse upwards towards 
the M point, where the surface states have an energy 
of about -0.8 eV. From M to X, this band disperses into 
the projected bulk band. The A^ band shows a small 
dispersion around M given rise to a large contribution on 
the Local Density of States (LDOS) of the first layer at 
an energy of about -1 eV, as shown in Fig. |31 

The A3 surface electronic states are at a lower energy 
than A^. The A3 states are due to the backbonds be- 
tween the anions localized in the first atomic layer and 
the cations in the second layer. The A3 band has a min- 
imum in the X high-symmetry point with an energy of 
-2.8 eV from the VBM. The band reaches its maximum 
at X' with an energy of -1.6 eV from the VBM. The band 
shows a dispersion of 1.2 eV, however, around X and M 
the band is almost flat, contributing to a large density 
of states in the first and the second layers at energies of 
about -2.8 eV and -2.5 eV, respectively. This can be ob- 
served on the LDOS in Fig.|31 where two peaks are found 
at these energies in the panels showing the LDOS in the 
first and the second layers. 

We found surface states with an energy of about 
—6.0 eV at the high-symmetry point X, that form a band 
denoted by C2. This band shows a large dispersion of 
about 2.5 eV, where the minimum of the band is at M 
with an energy of -6.3 eV, and its maximum is around 
X' with an energy of -3.8 eV. These surface states are lo- 
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cated at the cation (In) atoms and are due to the bonding 
between the In and As atoms at the top layer. From X to 
M, the band shows a small dispersion which is reflected 
in the LDOS where a large contribution is found at about 
-6.2 cV in the panel showing the projected LDOS in the 
first layer in Fig.|21 From M to X', the C2 band disperses 
upwards of about 2.5 eV given rise to a small contribution 
to the LDOS as shown in Fig. 13 

At lower energies we found another occupied surface 
electronic band denoted by A2 which extends almost 
along all the high-symmetry points in the 2DIBZ. These 
surface states are located in the anion (As) atoms, and 
have a s character due to the backbonds between the 
atoms at the first and second layers, and some contribu- 
tion is also found from the backbonds between the atoms 
in the second and third layers. From F to M going along 
X, the band does not show dispersion and is at -10.2 eV. 
From M to r going along X', the band has a dispersion 
of about 1 eV, showing a minimum around X'. 

We have also found several localized states at X be- 
tween r and M with an energy from -2 eV to -1 eV. These 
states are inside the projected bulk band, therefore, they 
are resonance-like states. These resonance states, de- 
noted by A4, disperse upwards from X towards both, 
F and M. Most of the A4 states have a p-character, 
and are localized at the anion in the first atomic layer. 
The A4 states with lower energy, between -1.9 eV and 
-2.0 eV, show also a s character and are localized at the 
third layer. 

Above the VBM wc found two unoccupied surface 
states bands, namely, C3 and C4. The C3 surface states 
are localized at the cations in the first and third atomic 
layers. They show a strong p character due to the dan- 
gling bonds at cations. At X, we found that C3 has a 
maximum with an energy of about 2 eV, and has its min- 
imum value between M and X' with a energy of about 
1.4 eV. Finally, at 2.7 eV from VBM we found empty 
surface states that form a band along all the high sym- 
metry points in the 2DIBZ. This band is denoted by C4, 
and show a very small dispersion along the 2DIBZ. 



IV. DISCUSSION 

In this section we discuss our results, and 
compare with available experimental measure- 
mentaSi9,i0iii,i2,i3,i4,i5.i6.i7,i8 ^nd theoretical calcu- 

lations^iiSiSii. The electronic properties of InAs(llO) 
have been investigated previously using experimental 
techniques like photoemission (pE)Si2iiS4iiii24i^ and in- 
verse photoemission (iPE)iiiiiiS4i& spectroscopies. We 
found no agreement between experimental measurements 
because they present difficulties to identify the position 
of the VBM or Ep, and the employed samples are quite 
different. Furthermore, photoemission measurements 
are difficult to interpret since emissions from surface 
states are usually hidden by emissions from bulk states 
in InAs surfaces, due to the small bulk-band gap. 



On the other hand, theoretical calculations have been 
performed from ab initio^ and semi-empirical TBi*^*S*i 
methods. We summarize our results and some of the 
experimental and theoretical data in Table ^ where 
we show the energy values at high symmetry points in 
2DIBZ of the surface electronic states denoted by A^, 
A3, A4 and C2. The first column shows our results, the 
next three columns show experimental measurements ob- 
tained by PB2i2iiS, and the last two columns show the- 
oretical results^i^. The values in Table ^ are those re- 
ported in the corresponding reference, or they have been 
estimated from the figures in each reference, then errors 
of about 0.1 eV in the estimated values are expected. 

Alves et alA performed an ab initio calculation based 
in the Density Functional Theory (DFT) within the Lo- 
cal Density Approximation (LDA), where many-body ef- 
fects were not taken into account. They considered slabs 
of only eight atomic layers (16 atoms), and the employed 
plane-wave basis set was expanded up to an energy cut- 
off of 8 and 18 Ry. They reported results for the equi- 
librium atomic structure and the electronic band struc- 
ture. It is known that the equilibrium atomic geome- 
tries can be found with good accuracy, but underestima- 
tion and/or overestimation of electronic states is always 
present in DFT-LDA calculations due to the approxi- 
mations employed, for example, usually DFT neglects 
many-body effects^^. Furthermore, plane-wave basis ex- 
pansions always present convergence problems to find lo- 
calized states. Therefore, the comparison of the surface 
states from ab initio calculations ^ with semi-empirical 
results and PE measurements always presents deviations 
up to ±1 eV. We also compare our results with semi- 
empirical TB calculations done by Mailhiot et alX. These 
TB calculations employed a theoretical method similar 
to the one used here but with different atomic positions, 
that were not fully relaxed. 

Both theoretical calculations^ found an empty surface 
state C3 and was identified with dangling bond states at 
cations. While the DFT calculation'^ found that C3 has a 
minimum at X, we obtained a maximum at the same sym- 
metry point in agreement with other semi-empirical cal- 
culations^i^iSii. DFT calculations reported an upwards 
dispersion from X of about 1.4 eV, while we found a 
downwards dispersion from X of about 0.6 eV. This dis- 
crepancy between DFT and semiempirical calculations is 
expected since DFT uses a small basis set that can not 
reproduce conduction states, while semiempirical calcu- 
lations with an extended basis set can do. The C3 surface 
states have been measured by using inverse photoemis- 
gjQj^i3ji4ji5ji6 pjjjy ^jjg high-symmetry point. 

The experimental measurements assigned an energy be- 
tween 1.7 eV and 1.9 eV at X, and evidence of an upwards 
dispersion from this point have been observecU^ in agree- 
ment with our calculations. However, a more detailed ex- 
perimental analysis is necessary to conclude more about 
empty states. 

The occupied surface states denoted by A5, A3, and 
C2, were also calculated in Refs. 3 and 4. Both calcula- 
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tions identified the surface states with the dangling 
bond states in the anion atoms, in agreement with our re- 
sults. The A5 surface state found using a first-principles 
method^t is shift 0.3 eV in average, above the experi- 
mental value^ and our calculation. Experimentally^, the 
A5 surface state has a dispersion of about 0.15 eV from 
X to X' through M, while we calculate a dispersion of 
0.5 eV, similar to the one obtained using DFTl. The 
surface states denoted by A4 have been also observed ex- 
perimentalljiSiS, and calculated by Mailhiot et alA. Like 
for the A^ states, we can not identify A4 states at the F 
point. PE measurements showed that these states have a 
dispersion of about 0.55 cV from X to to X' through M, 
which is in agreement with our calculated value of 0.5 eV, 
while Mailhiot et air found a smaller dispersion of 0.2 eV, 
and this value was not calculated using DFT'^. Previous 
semiempirical results'* also found the occupied states la- 
beled by A4, however, these surface states were not well 
identified since in their calculations the A4 states show a 
very similar dispersion than the A^ states. In general, the 
A4 are resonant states and they are difficult to calculate, 
specially if a plane-wave basis is used as the DFT calcu- 
lations discussed hereSi. Below A4, it has been observed 
other resonant states denoted by A3, using PEjSi. The re- 
ported data for A3 is quite different from our calculations 
and previous TB calculations^, maybe because the iden- 
tification of these states is not clear experimentally^'^. 
On the other hand, the C2 states can be experimentally 



identified since they are in a gap, except at the X point 
where they disperse into the projected bulk states. 



V. SUMMARY 

We performed a TB calculation using a fully relaxed 
atomic geometry to study the electronic structure of the 
clean InAs(llO) surface. A very detailed analysis was 
done, and a good agreement between our calculations 
and some experimental data was found. In conclusion, 
we found that fully relaxed atomic positions from ab 
initio methods in combination with our semi-empirical 
tight-binding calculation, better resembles photocmis- 
sion measurement of cleaveage InAs(llO) samples. Al- 
though the agreement between our results and other the- 
oretical and experimental data is good, we conclude that 
more experimental studies are necessary to clearly eluci- 
date the atomic relaxation and the electronic properties 
of InAs(llO). 



Acknowledgments 

We acknowledge the partial financial support from 
DGAPA-UNAM grant No. IN104201, CONACyT grants 
No. 36651-E and 36764-E. 



* Author to whom correspondence should be addressed. 
Emaircecilia@fisica.unam.mx" 

^ Group in Nitride Semiconductors Compounds: Physics 
and Applications, edited by B. Gil, (Oxford Science Pub- 
lications, Great Britain, 1998). 

^ See for example, F. Bechstedt and R. Enderlein, Semicon- 
ductor Surfaces and Interfaces, (Akademie-Verlag, Berlin, 
1988); A. Zangwill, Physics at Surfaces, (Cambridge Uni- 
versity Press, New York, 1992); H. Liith, Surfaces and In- 
terfaces of Solids, (Springer Verlag, Berlin, 1993) . 

^ J.L. Alves, J. Hebenstreit and M. Sheffler, Phys. Rev. B 
44, 6188 (1991). 

* C. Mailhiot, C. B. Duke and D. J. Chadi, Phys. Rev. B 
31, 2213 (1985). 

^ R. P. Beres, R. E. Allen and J. D. Dow, Phys. Rev. B 26, 
5207 (1982). 

^ B. Engels, P. Richard, K. Schroeder, S. Bliigel, Ph. Ebert 

and K. Urban, Phys. Rev. B 58, 7799 (1998). 
'' A. I. Shkrebtii, N. Esser, M. Kopp, P. Haier, W. Richter, 

and R. Del Sole, Appl. Surf. Sci. 104/105, 176 (1996). 
® C. B. M. Andersson, J. N. Andersen, P. E. S. Persson and 

U. O. Karlsson, Phys. Rev. B 47, 2427 (1992). 
® C. B. M. Andersson, J. N. Andersen, P. E. S. Persson and 

U. O. Karlsson, Surf. Sci. 398, 395 (1998) 
° D.M. Swantson, et al. Surf. Sci. 312, 361 (1994). 
^ D.M. Swantson, et al. Can. J. Phys 70, 1099 (1992). 
^ H. W. Richter, et al, J. Vac. Sci. Technol. B 4, 900 (1986). 
^ H. Carstensen, R. Claessen, R. Manzke, and M. Skibowski, 

Phys. Rev. B 41, 9880 (1990). 



^'^ W. Drube, D. Straub and F. J. Himpsel, Phys. Rev. B 35, 
5563 (1987). 

W. Gudat and D. E. Eastman, J. Vac Sci. Technol. 13, 831 
(1976). 

^® J. van Laar, A. Huiser and T. L. van Rooy, J. Vac Sci. 

Technol. 14, 894 (1977). 
" C. Mailhiot, C. B. Duke and D. J. Chadi, Surf. Sci. 149, 

366 (1985). 

^® C. B. Duke, A. Baton, A. Kahn and C. R. Bonapace, Phys. 
Rev. B 27, 6189 (1983). 

C. Noguez and S.E. UUoa, Phys. Rev. B 53, 13138 (1996); 
and references therein. 
^° P. Vogl, H. P. Hjalmarson and J. D. Dow, J. Phys. Chem. 
Solids 44, 365 (1983). 

T. J. Godin, J. P. LaFemina and C. B. Duke, J. Vac. Sci. 
Technol. A 10, 2059 (1992). 
^2 C. Noguez, Phys. Rev. B 58, 12641 (1998). 

"Electronic Structure and the Properties of Solids", Walter 
A. Harrison, (Dover Publications, Inc., New York, 1980). 
J. R. Chelikowsky and M. L. Cohen, Solid State Commun. 
29, 267 (1979). 

25 M.S. Hybertsen and S.G. Louie, Phys. Rev. B 34, 5390 
(1986); ihidZS, 4033 (1988); F. Bechstedt, and R. Del Sole, 
Phys. Rev. B 38, 7710 (1988); R. Godby, M. Schliiter and 
L.J. Sham, Phys. Rev. B 37, 10159 (1988). 



5 















ideal 


0.U4 


U.U 


6/4:00 


U.U 


do 


DPT* 


5.844 


0.70 


4.656 


0.122 


1.463 




5.861 


0.75 


4.663 


0.128 


1.445 


LEED 


6.036 


0.78 


4.985 


0.140 


1.497 






Lo (deg) 


(%)ciai 


(%) caai 


(%)cia2 


Ideal 


ao/2 


0.0 


0.0 


0.0 


0.0 


DPT* 


3.361 


30.7 


-1.80 


-0.22 


-2.00 




3.395 


32.0 


-1.18 


-0.18 


-1.82 


LEED 


3.597 


36.5 


-4.22 


+2.03 





TABLE L Structural parameters as defined in Fig.0 Param- 
eters obtained from DPT calculations using an energy cutoff 
of 8 Ry (*), and 18 Ry (**), both from Ref. 3, and parameters 
from Low-Energy Electron Diffraction (LEED) measurements 
from Ref. 12. 
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TABLE II: Experimental and theoretical values of the surface 
states at high-symmetry points of the 2DIBZ. The energy val- 
ues are in eV, where the zero energy corresponds to the VBM. 
(t) Estimated value from Pig. 7(a) in Ref. 3. (*) Estimated 
values from Pig. 9(b) in Ref. 4. An error of about 0.1 eV is 
expected in the estimated values. 



PIG. 1: Model of the atomic geometry of InAs(llO). (a) Top 
view of a surface unit cell, (b) Side view of the first three 
atomic layers of the surface, (c) Two-Dimensional Irreducible 
Brillouin Zone. 
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FIG. 2; Electronic band structure of the reconstructed 
InAs(llO) surface. Tiny dots represent the projected bulk 
states, while black dots represent surface electronic states. 



FIG. 3: Total and the projected local density of states in 
the first, second and third atomic layers of reconstructed 
InAs(llO). 
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